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Abstract 


In this report, we study the primal-dual algorithm as in [1] for 
the “uncapacitated facility location” problem. We explain the 3- 
approximation primal-dual approach for solving the problem and present 
the analysis of its worst-case approximation ratio. We also illustrate 
a sequence of asymptotically growing example sizes where the primal- 
dual approach approaches the approximation ratio of at most 3; we 
use the case of linear arrangements of two facilities, a singular client, 
and a cluster of clients that are almost equidistant from one of the 
facilities. We also prove that the approximation ratio is strictly less 
than 3, that is, it can never be equal to 3. For certain classes of in¬ 
puts, we illustrate that the nearest-neighbor approach works better, 
especially when for every facility i, f(i) < c[i][j], for all clients j. For 
some “trivial” cases, where facility costs for all the facilities are zero or 
negligible, we show that the primal-dual approach does not compute 
the optimal solution, whereas the nearest neighbor approach certainly 
computes the optimal solution in polynomial time. 

Next we introduced and explained a linear time algorithm (with re¬ 
spect to number of clients d) for a small constant number of facilities. 
The algorithm first finds by brute-force, all the possible combinations 
of facilities that can be opened and then using the nearest neighbour 
approach, assigns each client to its nearest opened facility. The combi¬ 
nation of facilities for which the total cost (facility cost of the opened 
facilities plus the service cost for all clients) is minimum, is the solu¬ 
tion. We further proved that the solution thus obtained is the optimal 
solution. Its time complexity is 0(d(f 2^ -1 )), which is 0(d) if / is a 
constant, that is linear with respect to demand size d. 

We implemented the primal-dual method, the nearest neighbour 
approach, and the linear time heuristic and performed certain exper¬ 
iments on classes of inputs where clients and facilities are randomly 
located in the plane (very rare chance of clustering of clients or fa¬ 
cilities). We compared the output of the primal dual approach and 
the nearest neighbour approach for better results. First we kept the 
facility cost much higher in comparison to the service cost c[i][j] and 
observed that for a small number of facilities and to a certain number 
of demands, the primal-dual approach gave better results, whereas 
even for a small value of /, when the demand d is comparatively 
higher, the nearest neighbour approach works better. For high val¬ 
ues of /, the nearest neighbour approach is better for high as well as, 
up to a certain point, lower values of d. Next we made facility cost 
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comparable to service cost, and observed that the nearest neighbour 
approach incurs less cost than primal dual method, always. 

Finally, we compared the time taken for computation of the out¬ 
put by the primal-dual method and the linear time heuristic. The 
motivation behind this experiment was to find a value of / for which 
computing the optimal solution is practically reasonable. With the 
help of the outputs we confirmed the linear dependence of running 
time on demand size for a constant value of / using the linear-time 
heuristic, and its exponential dependence on the number of facilities. 
We concluded that for a value less than 13 for /, it is computation¬ 
ally reasonable to use the linear-time approach over the primal-dual 
method for any size of demand and obtain the optimal solution. 


1 Introduction 

The uncapacitated facility location problem is a combinatorial optimization 
problem. It has a wide range of practical applications when it comes to the 
problems relating to setting up facility distribution centres. In this report we 
try to conclude how for different classes of this problem generated by changing 
the variables, a different set of approaches works better. This is important 
because not every time the cost of installing a facility is too huge, or the 
number of facilities that we want to install is high. Maybe we need to find 
better results if clients are clustered or when they are randomly distributed 
throughout the domain space. 


2 Survey of work in literature 

2.1 The uncapacitated facility location problem [1] 

In the uncapacitated facility location problem, we have a set of clients or 
demands D and set of facilities F. For each client j E D and facility i € F, 
there is a cost c tJ of assigning client j to facility i. Furthermore, there is a 
cost fi associated with each facility i £ F. The goal of this problem is to 
choose a subset of facilities F' C F so as to minimize the total cost of the 
facilities in F' and the cost of assigning each client j e D to the nearest 
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facility in F'. In other words, we wish to find F' C F such that, 

Minimize 22 & + 22 c v 

iEF' j£D min i£F' 

first part is called facility cost and the later one is called assignment cost or 
service cost. This is a computationally NP-hard problem until NP = P. 

Integer programming formulation for this problem. We define decision 
variables y* E {0,1} for each facility f t E F; if we decide to open facility 
i, then y t = 1 and y t = 0 otherwise. We also introduce decision variables 
Xij E {0,1} for all i E F and all j E D: if we assign client j to facility i, then 
Xij = 1 while = 0 otherwise. The objective function becomes, 

Minimize E fiVi + 22 c ij x ij 

iGF ifFjeD 

We need to make sure that each client j E D is assigned to exactly one 
facility. This can be done by, 

22 X ' J = i 

ieF 

We also need to make sure that the client is assigned to a facility that is 
open. This can be done by, 


Xij < Vi 

Thus, the integer linear programming formulation of the facility location 
problem can be summarized as follows: 


minimize 22 fiVi + 22 Ci i Xi i 


iGF ifFjeD 


subject to 22 Xi J = 

Vi e D, 

ieF 


VI 

Vi E F,j E D, 

Xij E {0,1}, 

Vi E F,j E D, 

Vi e {o,i}, 

i E F. 
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Linear programming relaxation from the integer linear programming can 
be obtained by replacing the constraint x t j E {0,1} and y, E {0,1} with 
Xij > 0 and y t > 0. Thus, the relaxed linear programming can be summarized 
as follows: 


minimize E fiVi + X] O, 

ieF if F,jeD 

subject to = 1, 

ieF 

J'ij E Ui i 
Xij E 
Vi > o, 


Vj e D, 

Vi E F,j E D, 
Vi E F,j E D, 
i E F. 


Forming dual of the relaxed linear programming: Ignoring the facility cost, 
i.e. to pretend that cost /* = 0, for all i E F, the optimal solution is to open 
all the facilities and assign each client to its nearest facility. We introduce 
a dual variable v 3 set to v 3 = min then this lower bound is v 3 . 

Now each facility takes its cost f) and shares it apportioned among the clients 
it provide service to: f t = J2jeD u ’y ’ where each Wi 3 > 0. A client j needs to 
pay this share only if it uses facility i. Hence, the optimal solution is to assign 
all the clients to the nearest facility, now we can set v 3 = min i(z F (ci 3 + w l3 ), 
and if we allow v 3 < c l3 + w l3 , the objective function maximizing ]W e£) 'en¬ 
forces Vj to be equal to the smallest right-hand side overall facilities i E F, 
then any feasible solution to dual is lower bound on cost of optimal value on 
facility location problem. Thus, the dual linear program for the primal linear 
program can be summarized as: 

maximize v 3 
jeD 

subject to E Wij < fi, Vi, E F 

jeD 

Vj - w i3 < Ci 3 , Vi E F,j E, D 

>0, Vi E F,j E D. 
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2.1.1 3-approximation primal-dual algorithm for the uncapaci¬ 
tated facility location problem [1] 

Let (v*, w *) be a feasible solution to dual which is maximal, that is v* cannot 
be increased by any positive amount and then derive a set of wf that gives 
a dual feasible solution. 

Definition 1. [1] Given a dual solution (v*,w*), we say that a client j 
neighbours facility i if v* > c i3 . We let N(j) = [i £ F : v* > c i3 } be the 
neighbours of a client j and N (i) — {j £ D : v* > Ci 3 } be the neighbours of 
a facility i. 

Definition 2. [1] Given a dual solution (v*,w*), we say that a client j 
contributes to a facility i if tvf > 0. 

The primal-dual algorithm is as follows: 

The idea is to generate a maximal dual solution by increasing the dual vari¬ 
able Vj. Let S be the set clients whose dual variables are being increased 
and let T be the set of facilities whose dual inequality is tight. Initially, set 
S = D and T = <f. increase v 3 uniformly for all j £ S. Once v 3 = c l3 for some 
i, increase Wij uniformly with Vj until one of the two things happen: either j 
becomes neighbour of some facility in T, or a dual inequality becomes tight 
for some facility i. In first case, remove j from S', and in second case add i 
to T, and remove all the neighbouring clients N(i) from S. Once S is empty 
and every client neighbours a facility in T, select a subset T' of T such that 
no client is neighbour to more than one facility in T'. This can be done 
by iterative picking of an arbitrary facility i £ T, and adding it to T' and 
then deleting all the facilities i' £ T such that there exist some client j that 
contributes to both i and i!. once its done then open all facilities in V and 
assign every client to the closest open facility. 

Lemma 1. [1] If a client j does not have a neighbour in T', then there exist 
a facility i £ T' such that c^j < 3 Vj. 


Proof: Let j be an arbitrary client that does not neighbour a facility in 
T'. During the course of the algorithm, v 3 stopped increasing because j 
neighboured some h £ T, obviously h f T', since otherwise j would neighbour 
a facility in T'. the facility h must have been removed from T' during the 
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Figure 1: Proof for lemma 1 

final subset picking because there exist some other client k that contributes 
to both h and some facility i e T'. Now, Chj < v 3 because j neighbours h. 
At the point in algorithm at which Vj stopped increasing, by the choice of 
h, at this point either h is already in T or algorithm adds h to T. Now that 
client k contributes to facility h, it must be the case that either ty has already 
stopped increasing or at same point that v 3 stopped increasing. Because the 
dual variables are also increase uniformly, we have v 3 > ty. As k contributes 
to bot, facility i and h, Vk > Chk and v + k > c^. Thus, Vj > Vk > Chk and 
U j ^ ^ k ^ Cik • 

Using triangle inequality c l3 < Cik + Ckh + Chj• This implies that c l3 < 3 v 3 

Theorem 1. [1] Above mentioned algorithm is a 3-approximation algorithm 
for uncapacitated facility location problem. 

Proof: For ay client that contributes to a facility in T', assign it to this 
facility. Note that by construction of the algorithm, any client contributes 
to at most one client in T 1 , so this assignment is unique, for clients that 
neighbours facilities in V but do not contribute to any of them, assign each 
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some arbitrary facility in T'. Let A(i) C N(i) be the neighbouring clients 
assigned ti the facility i E T'. Then as discussed above the cost of opening 
the facilities in V plus the cost of assigning the neighbouring clients is, 

+ E C; ) = E E ("o + °ij) = 12 12 r b 

i£T' jEA(i) i£T' j£A(i) i£T / j£A(i) 

where the first equality holds because i G V implies Wij = f, and 

Wij > 0 implies.; in A (7). Let Z be the set of all clients not neighbouring a 
facility in T', so that Z = D — Uj G T'A(i). We have by Lemma 1 that the 
cost of assigning j e Z to some facility in V is at most 3 Vj. Thus the total 
assignment cost for these clients is at most 3 J2jez £■ 

Example: 


2.2 Best results in Literature 

Following are the works done in this field. 

3.16 - approximation algorithm by D. Shmoys, E. Tardos and K. Aardal [2], 
They used filtering and rounding algorithms build on the results of Lin and 
Vitter [7] [8] who gave an elegant technique of filtering for fractional solutions 
to linear programming relaxations and used it to give O ( logn )-approxiination 
algorithm for uncapacitated facility location problem. And showed that for 
symmetric assignment costs in the k-median problem following triangle in¬ 
equality, one can find a solution cost of 2 (i + e) times the optimum. 

2.408 - approximation algorithm by S. Guha and S. Khuller [3]. They com¬ 
bined a greedy heuristic with the algorithm by D. Shmoys, E. Tardos and 
K. Ardal [2] and also proved that there can be no better than 1.463 approx¬ 
imation ratio for metric uncapacitated facility location problem using the 
reduction from set cover problem. 

1.736 - approximation algorithm by F. Chudak [4], They obtained a (1+2/e)- 
approximation algorithm where 1 + 2/e ~ 1.736, by rounding an optimal 
fractional solution to a linear programming relaxation and generalization of 
decomposition techniques of Shmoys, Tardos and Aardal [2], 
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1.582 - approximation algorithm by Maxim Sviridenko [5]. they used clus¬ 
tering techniques and rounding technique for the problems with cardinality 
constraints. They used theses techniques to derandomize the randomized 
algorithms without actually constructing them. 

1.52 - approximation algorithm by M. Mahdian, Y. Ye and J. Zhang [6]. 
They used the greedy algorithm of Jain, Mahdian and Saberi [9] combined 
with the greedy augmentation of Charikar, Guha and Khuller [3] [10], and 
showed that uncapacitated facility location problem can be approximated 
within a factor of 1.52. 


3 Example that is approximately 3-OPT 





U_1 

u_2 

u_i 
u n 


Figure 2: Facility Location Example 


In the example shown in figure 2 , there are two facilities /i and f 2 and 
n + 1 clients v ly iq, u 2 ,...., u n where the n clients iq, u 2 , ...u n are clustered 
equidistant from facility / 2 and almost equidistant from /i. 

The facility costs of /i = 1 and f 2 = 3. The c tJ cost of the following facility 
client arrangement is as follows. 


Vi Ui 

/i k 3k + 1 

f 2 k k + 1 


u 2 u 3 

3k+ 1 3k+ 1 
k + 1 k + 1 


3k + 1 
k + 1 
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Algorithm will increase the Vj for all the clients to 5 and hits minimum c V] 
with client v\ and it becomes neighbour to both facilities f\ and / 2 . As the 
set T is empty, v\ cannot become neighbour to any client in T. hence while 
gradually increasing the Vj and w t] , inequality becomes tight for fi and it is 
put in T and are removed from S. See figure 3. In next iteration, Vj are 
raised to 6 as Vj = c Ul2 = 6. All u, become neighbour to / 2 . Now gradually 
increasing the Vj and w l; j , inequality becomes tight for / 2 before u\ becomes 
neighbour to f\. fo is added to T and all iq are removed from S. See figure 
4. 



Figure 3: Facility Location Example 



Figure 4: Facility Location Example 


As all clients in S have become neighbour to at least one facility in T, the 
algorithm now picks T'. It starts picking first the facility that was selected 
first in T and picked /i to begin with and puts it in T'. Now client v\ is 
connected to both f\ and / 2 , hence / 2 cannot be added to T', therefore all 
the clients that are neighbour to / 2 are connected to /i. see figure 5. 
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This facility location assignment costs total of (1 + k + (3 k + l)n), see 
figure 6, whereas the Optimal solution is to join {«i, m 2 , u n } to /2 and 
v\ to /i, costing 4 + k + (k + 1 )n, see figure 7. As k and n increase, the 
approximation ratio tends to converge to 3. 


fl 


V_1_ 


U_1 

u_2 

u _i 
u n 


Figure 5: Facility Location Example 


fl 


v 1 




U_1 

u_2 

u _i 
u n 


Figure 6: Facility Location Example 



Figure 7: Facility Location Example (OPT) 
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4 Primal-dual vs nearest neighbour approach 

4.1 For facility costs equals to Zero 

Nearest neighbour algorithm is bound to give the optimal solution in 
this case, because the cost of installation of a facility is Zero, the optimal 
solution is to open all the facilities and connect every client to the nearest 
facility. 

Primal-Dual algorithm may or may not give the optimal solution. That 
is so because while building set T it allows a client to be neighbour to more 
than one facility and then while picking up the clients for T' it may close 
some facility due to some conflicting client, hence increasing the cost of as¬ 
signment of clients to facilities. 

Example:: 


fl v_l f2 v 2 



Figure 8: Example with zero facility cost 

Consider figure 8, there are two facilities and two clients, both the facility 
has facility establishment cost Zero, and following c l3 values: 


Vi v 2 

/i k 2k+ 1 

f 2 k k + 1 


As the Vj increases, it becomes equal to k for client v\ with facilities f\ and 
/ 2 both. V\ becomes neighbour to both the facilities. Because of Zero facility 
cost, inequality becomes tight for both the facilities but chosen randomly, /i 
is selected and added in T, and v\ is removed form S. Vj are again increased 
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to k + 1 and it becomes equal to c[/ 2 ][u 2 ], inequality for / 2 becomes tight 
before u 2 could neighbour fi. hence / 2 is chosen and added to T and u 2 is 
removed from S. 


fl v_l f2 

■ - + -■ 


v 2 



Figure 9: Example with zero facility cost 


fl 


v_l f2 

- ■ 


6 


Figure 10: Example with zero facility cost 

While selecting facilities in T' , /j is chosen first, now as V\ is neighbour 
to both fi and / 2 , / 2 cannot be chosen as a facility in T'. because a client 
can neighbour to only one facility, hence client u 2 is connected to f\ costing 
2k + 1. 

fl - >d -- f2 v _2 

on o 


Figure 11: Example with zero facility cost 

Cost of providing service using nearest-neighbour approach would have 
been 2k + 1, by joining v\ to one of the any facility and u 2 to / 2 . Whereas, 
using primal-dual approach , it costs total of 3k + 1 which is more than OPT 
= 2k + 1. 
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4.2 Linear time complexity for constant number of fa¬ 
cilities. 

Services like petrol pumps stations are generally in a very limited number in 
a region, which also remains constant for a long time and only the number 
and arrangement of clients keep on changing. In this section we are trying 
to find an algorithm, linearly complex with respect to d, that is the number 
of clients with a fixed number of facilities, for the facility location problem 
using the approach of nearest neighbour. 

4.2.1 / = 3 

In this case, the number of facilities is 3, i.e. / = 3. We form all possible 
combinations of facilities, T\ for the k th combination. There will be 7 such 
combinations: 

Case 1: 1 combination when all three facilities are chosen. 

Case 2: 3 combinations when two out of three facilities have been chosen. 
Case 3: 3 combinations when one out of three facilities have been chosen. 

In case 1, for every client i, join it with the facility j such that c %J is 
minimum over all facilities. 

In case 2, for each combination, for every client i, join it with the facility 
j such that Cij is minimum over the two facilities. 

In case 3, for each combination, join every client i with the facility j. 

For all the 7 combinations, compute the total cost, that is the sum of the 
facility cost of all the facility in set T\ and the c tJ values of the facility-client 
connection for all the clients. 

The combination set T[ with minimum total cost is the solution. 

Analysis For every combination, every client is picking the minimum c l7 
which can be, using simple if statement, done in 0(1) time. Hence there will 
be n such computations, one for each client. Hence, every combination is of 
0(n) complexity. And there are 7 such combinations. Therefore the total 
complexity of the algorithm is 7 * 0(d), that is 0(d). 
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4.2.2 / = 4 


In this case, the number of facilities is 4, i.e. / = 4. We form all possible 
combinations of facilities, T\ for the k th combination. There are 15 such 
combinations: 

Case 1: 1 combination when all four facilities are chosen. 

Case 2: 4 combinations when three out of four facilities have been chosen. 
Case 3: 6 combinations when two out of four facilities have been chosen. 
Case 4: 4 combinations when one out of four facilities have been chosen. 

In case 1, for every client i, connect it to the facility j such that c V) is the 
minimum over all facilities. 

In case 2, for each combination, for every client i, connect it to the facility 
j such that Cij is the minimum over the the three facilities. 

In case 3, for each combination, for every client i, connect it with the 
facility j such that c tJ is the minimum over the two facilities. 

In case 4, for each combination, connect every client i with the facility j. 

For all the 15 combinations, we compute the total cost, that is the sum 
of the facility cost of all the facilities in set T\ and the c tJ values of the 
facility-client connection for all the clients. 

The combination set T' k with minimum total cost is the solution. 

Analysis For every combination, every client is picking the minimum c tJ 
which can be, using simple if statement, done in 0(1) time. Hence there will 
be n such computations, one for each client. Hence, every combination is of 
0(d) complexity. And there are 7 such combinations. Therefore the total 
complexity of the algorithm is 15 * 0(d), that is 0(d). 

4.2.3 Extending to general case. 

Let / be the number of facilities and n the number of clients. We form all 
possible combinations of facilities, T\ for the k th combination. 

Case x: i is the number of facilities being picked. There will be ?C X such 
combinations. For every combination in this case, every client will need x 
comparisons to find the minimum c l3 . 
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Computations required in case x =* C x xd. 

f 

Total computations required = £ fC x xd. 

X=1 

f 

= dY, f C x x. ...( 1 ) 


Using Binomial expansion... 

(x + l) n = n C 0 x°r + n Cxx 1 !™- 1 + n C 2 x 2 l n ~ 2 + ... + n C n x n l° 
Differentiating both sides w.r.t. x 

n(x + If 1 = 0 + 1 * n Cix°l n ~ l + 2* n C 2 x l l n ~ 2 + ... + n * n C n x n ~ l l° 
Substituting x — 1 

n2 n ~ 1 = 1 C x + 2 * n C 2 + 3 * n C 3 + 4 * n C 4 + ... + n * n C n 

n 

n2 n ~ 1 = ^2 n C x x ...( 2 ) 

X=1 

from eq. (1) and (2)... 

Total computation required = d(f 2^ -1 ) 

Time complexity = 0(d(f2^~ 1 )) 

Out of the 2-f — 1 combinations, one with the minimum cost is the solution. 

If the number of facilities / is a constant then the time complexity becomes 
0(d). In general, the time complexity is 0(df 2-f -1 ), which is exponential in 
/• 

4.2.4 Proof that the solution obtained is optimal. 

Let in the optimal solution OPT set T of facilities is opened and every client 
is connected to exactly one of the facilities in T. We can say that exactly 
one of the combination sets T( is equal to T let that be T(,°. And that is so 
because T' k are all the possible combinations of facilities and T has to be one 
of them. While joining clients to facility in the T[°, every client i is joined 
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to the facility j l nearest to it. Suppose in OPT for some client i° this is not 
true and it is joined to some other facility j / j\ Now c ( o y > c i0 ji , by the 
dehnition of c l0] i. If it is equal then the algorithms solution is also a valid 
optimal solution, else cost of algorithm’s solution will be less than that of 
OPT which is not possible, by the definition of OPT. Hence the facility 
client connection for set T^° is optimal. 

Even for any other combination total cost cannot be smaller than that of 
OPT or else that combination will be OPT. 
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5 Implementation and Experimentation 

5.1 Performance of the primal-dual scheme Vs nearest- 
neighbour approach 

To estimate the approximation ratio for the runs, we need to estimate a 
lower bound for the optimal measure OPT, knowing that it is NP-hard to 
compute and optimal solution with measure OPT. So, instead we used a 
trivial lower bound on OPT rather than OPT itself, let C, be the minimum 
Cij value for each client , i.e. for i th client at least C t cost is added in the sum. 
let f° be the minimum facility cost among all. Now any optimal solution has 
to open at least one facility, and has to join every client to exactly one facility. 
Hence, 


OPT>f°pJ2°i 

Being a lower bound to optimal, the approximation ratio computed can 
be greater than the actual ratio. 

Note : Number of facilities (/) and number of clients (d) are taken as user 
input and all the other necessary inputs are randomly generated by the code 
within the mentioned bounds. The service cost (c[i][j]) is computed as the 
euclidean distance between client j and facility /. 


5.1.1 Facility cost is much higher than C[i] [j] 

Facility cost G [20 d, 21 d] 

X coordinate of client G [0, d] 

Y coordinate of client G [0, d] 

X coordinate of facility G [0, d] 

Y coordinate of facility G [0, d] 

where d is the number of clients. 
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Result Table 1 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

3 

10 

230 

246 

1.069 

438 

3 

20 

508 

590 

1.161 

1340 

3 

30 

969 

1091 

1.125 

2201 

3 

50 

1796 

2332 

1.298 

3826 

3 

100 

6806 

8681 

1.275 

10880 

3 

150 

11611 

12522 

1.078 

17815 

3 

200 

18359 

21172 

1.153 

26699 

3 

400 

64535 

100976 

1.564 

80909 

3 

750 

260445 

400888 

1.539 

291129 


Total Cost vs Demand (F=3) 


2.0DDD 



■ lower bound 


■ primal dual 


— •—nearest neighbour 


Figure 12: Total cost vs Demand (F = 3) 

For the first seven outputs, the total cost computed by the primal-dual 
method is smaller than the nearest-neighbour approach (figure 12). Whereas 
as the client size d increases, the cost of primal-dual method becomes higher 
than that of nearest-neighbour approach (figure 13). 
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Total Cost vs Demand (F=3) 
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Figure 13: Total cost vs Demand (F = 3) 


2 

1.5 

1 

0.5 

0 
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Figure 14: Approximation ratio (f = 3) 

This is so because, even when the number of clients were small, irre¬ 
spective to the cost of facility, nearest neighbour approach would open more 
facilities to provide service to nearest clients. Whereas primal-dual method 
opens less number of facilities and join far away clients with already opened 
facilities. But, as the number of clients increases, the domain space of loca¬ 
tion for clients also increases, forcing the primal-dual method to open more 
facilities hence increasing the cost. 


ratio (F = 3) 
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Result Table 2 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

10 

10 

219 

240 

1.095 

1045 

10 

25 

590 

778 

1.318 

4696 

10 

50 

1471 

2386 

1.622 

8666 

10 

100 

3662 

7980 

2.179 

22221 

10 

250 

15633 

34606 

2.213 

61628 

10 

500 

49309 

138195 

2.802 

141999 

10 

1000 

190262 

493006 

2.591 

376567 

10 

2000 

706677 

1997930 

2.827 

1079240 
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Figure 15: Total cost vs Demand (F = 10) 
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Total Cost vs Demand (F=10) 



* lower bound * primal dual • nearestneighbour 


Figure 16: Total cost vs Demand (F = 10) 


ratio (F = 10) 



Figure 17: Approximation ratio (f = 10) 
Trend similar to the previous one can be seen here. 
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Result Table 3 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

25 

25 

571 

732 

1.281 

7729 

25 

50 

1252 

2187 

1.746 

18798 

25 

100 

2929 

6413 

2.189 

49921 

25 

200 

7967 

24058 

3.019 

106879 

25 

500 

34196 

116702 

3.412 

281375 

25 

750 

79571 

328494 

4.128 

449530 

25 

1500 

296401 

1383900 

4.669 

1034000 

25 

2500 

821446 

3397940 

4.136 

2052550 

25 

4000 

2132010 

7767130 

3.643 

4111410 


500000 
400000 
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■g 200000 

I— 

100000 
0 

0 


Total Cost vs Demand (F=25} 
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►nearest neighbour 


Figure 18: Total cost vs Demand (F = 25) 

Once again, the cost by primal-dual method and nearest neighbour follows 
the same trend where primal-dual cost is lower for less number of clients and 
higher than nearest neighbour approach for higher number or clients. 
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Total Cost vs Demand (F=25) 

10000000 
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-primal dual 


- nearest neigh hour 


Figure 19: Total cost vs Demand (F = 25) 


ratio (F = 25) 

5 



0 1 — 

0 1000 2000 3000 4000 5000 

Figure 20: Approximation ratio (f = 25) 

Point to be noticed here is that even though we are using 3-approximation 
algorithm, the approximation ratio in figure 20 exceeds 3. This is so because 
we have used a very loose lower bound to optimal for calculating the ratio, 
and as mentioned in the beginning of this section, the ratio may exceed 3. 
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Result Table 4 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

50 

50 

1184 

1923 

1.624 

29814 

50 

100 

2703 

8564 

3.168 

88607 

50 

250 

9459 

39676 

4.194 

244792 

50 

500 

29571 

193034 

6.527 

520556 

50 

1000 

96566 

541525 

5.607 

1103500 

50 

2000 

364922 

2054000 

5.628 

2370500 

50 

3000 

698625 

4413940 

6.318 

3710200 

50 

5000 

1909600 

12924200 

6.768 

6932110 

50 

7500 

4267470 

29858600 

6.996 

11774500 

50 

10000 

7832380 

57699300 

7.366 

17896400 
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Figure 21: Total cost vs Demand (F = 50) 
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Total Cost vs Demand (F=50) 



-lower bound 
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Figure 22: Total cost vs Demand (F = 50) 


ratio (F = 50) 



Figure 23: Approximation ratio (f = 50) 
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Result Table 5 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

80 

80 

2031 

4266 

2.100 

66077 

80 

200 

6366 

29014 

4.557 

277033 

80 

500 

23824 

141815 

5.952 

814231 

80 

1000 

79323 

564541 

7.116 

1682110 

80 

2000 

269476 

2607100 

9.674 

3517810 

80 

4000 

983564 

9787310 

9.950 

7449750 

80 

7500 

3536770 

24302500 

6.871 

15687800 

80 

12000 

8377380 

55677400 

6.646 

27819600 

80 

20000 

24750700 

232214000 

9.382 

57156300 
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2000000 

0 
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Figure 24: Total cost vs Demand (F = 80) 
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Total Cost vs Demand (F=80) 



lower bound ■ primal dual ■ nearest neighbour 


Figure 25: Total cost vs Demand (F = 80) 


ratio (F = 80) 



Figure 26: Approximation ratio (f = 80) 

We see that the value of approximation ratio is increasing with the number 
of facilities. For / = 25 maximum ratio was 4.669, for / = 50 maximum ratio 
was 7.366, for / = 80 maximum ratio is 9.95 and so on. 


28 








Result Table 6 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

125 

125 

3141 

9588 

3.052 

190311 

125 

250 

7787 

30561 

3.924 

499513 

125 

500 

21231 

128072 

6.032 

1160580 

125 

1000 

66142 

435838 

6.589 

2567230 

125 

2000 

232640 

1917500 

8.242 

5242250 

125 

5000 

1385100 

17740000 

12.80 

14088400 

125 

9000 

3942840 

34806700 

8.827 

26787600 

125 

15000 

10822000 

89425200 

8.263 

49049900 


Total Cost vs Demand (F=125) 
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Figure 27: Total cost vs Demand (F = 125) 


29 



















Total Cost vs Demand (F=125) 
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Figure 28: Total cost vs Demand (F = 125) 


ratio (F = 125) 



Figure 29: Approximation ratio (f = 125) 

This continuous increase in the approximation ratio with the increasing 
number of facilities can be explained by the fact that the lower bound adds 
only one facility cost that too the one which is minimum. As the number of 
facilities increases and the client size increases, the number of facilities to be 
opened in the optimal solution is far more than just one. And as the facility 
costs are high, the increase in the approximation ratio is also significantly 
higher with the increase in number of facilities. 
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Result Table 7 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

250 

250 

7053 

30036 

4.258 

729808 

250 

500 

18155 

124654 

6.866 

2059710 

250 

1000 

52193 

562765 

10.78 

4725350 

250 

1500 

101774 

1086910 

10.67 

7627920 

250 

2500 

255370 

4445100 

17.40 

12869400 

250 

4000 

594090 

8324370 

14.01 

21008200 

250 

6500 

1463480 

24288100 

16.59 

34626200 

250 

10000 

3434700 

59388400 

17.29 

54470600 

250 

17000 

9503970 

124226000 

13.07 

196200400 
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Figure 30: Total cost vs Demand (F = 250) 
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Figure 31: Approximation ratio (f = 250) 
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Result Table 8 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

500 

500 

15588 

132653 

8 

3088100 

500 

2500 

191206 

3349320 

17 

24559100 

500 

1000 

44124 

599754 

13 

8091230 

500 

5000 

676103 

12668600 

18 

51117600 

500 

6500 

1074890 

24786800 

23 

67428900 

500 

9000 

2089440 

43020200 

20 

93959200 

500 

12000 

3668220 

66148600 

18 

12638800 


Total Cost vs Demand (F=500) 
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Figure 32: Total cost vs Demand (F = 500) 
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ratio {F = 500) 
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0 2000 4000 6000 8000 10000 12000 14000 

Figure 33: Approximation ratio (f = 500) 

For / = 250 and / = 500 primal dual approach gives a better solution 
for small as well as large set of clients than the nearest neighbour approach. 
Even though number of facilities is large, the probability of opening many 
facilities by the primal-dual approach is low, whereas almost all facilities 
opened by nearest neighbour approach. Even though the nearest neighbour 
approach uses minimum c[i][j\ value for each client, the high value of facility 
cost increases its total cost in comparison to primal-dual method. 
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5.1.2 Facility cost is similar to C[i] [j] 

Facility cost G [d, 2d] 

X coordinate of client G [0, 3d] 

Y coordinate of client G [0, 3d] 

X coordinate of facility G [0, 3d] 

Y coordinate of facility G [0, 3d] 
where d is the number of clients. 

In all the output sets for different number of facilities tabulated below, 
the total cost incurred by nearest neighbour approach is smaller than the cost 
incurred by the primal-dual approach. This is so because the facility cost for 
this series of experimentation is kept low that is very much comparable to the 
c [i][j] values. Keeping low facility cost forces even the primal-dual approach 
to open many facilities even if the input sizes are small. Moreover, as the 
service cost (c[i][j]) are similar to facility cost, indirectly joining a far away 
client to an opened facility proves to be much more expensive than opening 
the closest facility and connecting clients to it. 

The trend of approximation ratio, increasing with the number of facilities 
is still the same as before. Reason being the same; loose lower bound on 
OPT is being used to calculate the ratio. 
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Result Tables 1 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

3 

10 

88 

127 

1.443 

144 

3 

20 

245 

471 

1.922 

361 

3 

30 

659 

916 

1.389 

815 

3 

50 

1654 

2470 

1.493 

1925 

3 

100 

6102 

11992 

1.965 

6673 

3 

150 

14202 

17493 

1.231 

14988 

3 

200 

21243 

46274 

2.178 

22379 

3 

400 

86777 

131715 

1.517 

88848 

3 

750 

287800 

512893 

1.782 

291701 


Total Cost vs Demand (F = 3) 
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Figure 34: Total cost vs Demand (F = 3) 
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Figure 35: Total cost vs Demand (F = 3) 


ratio (F = 3) 



Figure 36: Approximation ratio (f = 3) 
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Result Tables 2 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

10 

10 

48 

93 

1.937 

141 

10 

25 

300 

498 

1.66 

716 

10 

50 

867 

2738 

3.158 

2002 

10 

100 

3612 

9432 

2.611 

5663 

10 

250 

24682 

63388 

2.568 

29917 

10 

500 

72178 

213702 

2.960 

82913 

10 

1000 

341155 

794783 

2.329 

361876 

10 

2000 

1.21E+06 

2.46E+06 

2.028 

1.25E+06 


Total Cost vs Demand (F = 10) 
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Figure 37: Total cost vs Demand (F = 10) 


















Total Cost vs Demand {F = 10} 
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Figure 38: Total cost vs Demand (F = 10) 


ratio (F = 10) 



Figure 39: Approximation ratio (f = 10) 
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Result Tables 3 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

25 

25 

195 

624 

3.2 

1008 

25 

50 

615 

2013 

3.273 

3224 

25 

100 

2231 

8858 

3.970 

7985 

25 

200 

8480 

36384 

4.290 

20712 

25 

500 

53212 

231398 

4.348 

82002 

25 

750 

117139 

439689 

3.753 

161130 

25 

1500 

485210 

1.27E+06 

2.622 

574860 

25 

2500 

1.37E+06 

4.38E+06 

3.198 

1.51E+06 

25 

4000 

3.39E+06 

1.20E+07 

3.553 

3.63E+06 


Total Cost vs Demand (F = 25) 
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Figure 40: Total cost vs Demand (F = 25) 
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Total Cost vs Demand (F = 25} 
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Figure 41: Total cost vs Demand (F = 25) 


ratio (F = 25) 



0 500 1000 1500 2000 2500 3000 3500 4000 4500 


Figure 42: Approximation ratio (f = 25) 


4500 
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Result Tables 4 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

50 

50 

458 

3355 

7.325 

3480 

50 

100 

1675 

11683 

6.974 

9850 

50 

250 

9873 

52876 

5.355 

38568 

50 

500 

40962 

264230 

6.450 

100982 

50 

1000 

145540 

691617 

4.752 

269056 

50 

2000 

593558 

2.91E+06 

4.898 

838253 

50 

3000 

1.34E+06 

5.37E+06 

4.012 

1.71E+06 

50 

5000 

3.92E+06 

1.94E+07 

4.950 

4.55E+06 

50 

7500 

8.66E+06 

3.54E+07 

4.089 

9.57E+06 

50 

10000 

1.44E+07 

5.66E+07 

3.941 

1.56E+07 


Total Cost vs Demand {F = 50) 
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Figure 43: Total cost vs Demand (F = 50) 
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Figure 44: Total cost vs Demand (F = 50) 


ratio (F = 50) 
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Figure 45: Approximation ratio (f = 50) 
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Result Tables 5 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

80 

80 

956 

5900 

6.171 

9895 

80 

200 

4870 

43590 

8.950 

39488 

80 

500 

29299 

217342 

7.418 

125858 

80 

1000 

114240 

817788 

7.158 

308772 

80 

2000 

442259 

3.45E+06 

7.810 

834688 

80 

4000 

1.87E+06 

1.24E+07 

6.655 

2.67E+06 

80 

7500 

6.14E+06 

3.02E+07 

4.916 

7.62E+06 

80 

12000 

1.66E+07 

8.42E+07 

5.058 

1.90E+07 

80 

20000 

4.72E+07 

1.97E+08 

4.179 

5.11E+07 


Total Cost vs Demand {F = 80) 
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Figure 46: Total cost vs Demand (F = 80) 
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Total Cost vs Demand (F = 80) 
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Figure 47: Total cost vs Demand (F = 80) 
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Figure 48: Approximation ratio (f = 80) 


ratio (F = 80) 



45 







Result Tables 6 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

125 

125 

1703 

15825 

9.292 

22452 

125 

250 

6194 

60262 

9.729 

67781 

125 

500 

25259 

248834 

9.851 

165761 

125 

125 

1665 

19513 

11.71 

24514 

125 

250 

6330 

62458 

9.866 

69006 

125 

500 

23937 

201810 

8.430 

163030 

125 

1000 

100155 

785474 

7.842 

404172 

125 

2000 

377345 

2.97E+06 

7.859 

998568 

125 

5000 

2.21E+06 

1.81E+07 

8.186 

3.75E+06 

125 

9000 

7.62E+06 

4.56E+07 

5.980 

1.04E+07 

125 

15000 

2.16E+07 

1.24E+08 

5.741 

2.62E+07 
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Figure 49: Total cost vs Demand (F = 125) 
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Figure 50: Total cost vs Demand (F = 125) 
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Figure 51: Approximation ratio (f = 125) 
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Result Tables 7 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

250 

250 

4639 

51904 

11.18 

93800 

250 

500 

18488 

298925 

16.16 

260728 

250 

1000 

66506 

1.01E+06 

15.25 

658120 

250 

1500 

152803 

1.71E+06 

11.21 

1.04E+06 

250 

2500 

402398 

5.56E+06 

13.80 

1.95E+06 

250 

4000 

1.05E+06 

1.05E+07 

10.07 

3.51E+06 

250 

6500 

2.72E+06 

2.87E+07 

10.57 

6.75E+06 

250 

10000 

6.67E+06 

4.85E+07 

7.276 

1.28E+07 

250 

17000 

1.85E+07 

1.35E+08 

7.338 

2.90E+07 


Total Cost vs Demand (F = 250) 
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Figure 52: Total cost vs Demand (F = 250) 
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Figure 53: Approximation ratio (f = 250) 


49 





















Result Tables 8 


f 

d 

lower 

bound 

Primal- 

dual 

ratio 

nearest 

neighbour 

500 

500 

12484 

204641 

16.39 

372261 

500 

1000 

50576 

1.01E+06 

19.88 

1.03E+06 

500 

2500 

289489 

5.76E+06 

19.91 

3.19E+06 

500 

5000 

1.19E+06 

2.50E+07 

21.09 

7.24E+06 

500 

6500 

1.97E+06 

3.12E+07 

15.83 

9.96E+06 

500 

9000 

3.72E+06 

5.25E+07 

14.12 

1.50E+07 

500 

12000 

6.54E+06 

1.01E+08 

15.49 

2.14E+07 
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Figure 54: Total cost vs Demand (F = 500) 
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Figure 55: Approximation ratio (f = 500) 
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5.2 Primal-dual scheme vs Linear-time heuristic 

In this section we analyze and compare the running time for the primal- 
dual scheme and the heuristic proposed in section 4.2 which is linear in d 
by simulating the codes for both schemes for different number of facilities 
and clients. We know that the primal-dual approach is an approximation 
algorithm where as the linear time heuristic gives the optimal solution. 

The motive behind this comparison is to find the upper bound in the 
number of facilities up to which it is computationally beneficial to use linear 
time heuristic because it is meant for a small constant number of facilities /. 
It is linear only with respect to number of clients d, computation increases 
exponentially with the increase in number of facilities for linear time heuristic. 

Number of facilities / and number of clients/demands d are taken as user 
input and the rest of the inputs required are randomly generated by the code. 
Service costs C[i\ [j] are computed as the euclidean distance between client j 
and facility i. 

Facility cost £ [d, 2d] 

X coordinate of client £ [0, 3d] 

Y coordinate of client £ [0, 3d] 

X coordinate of facility £ [0, 3d] 

Y coordinate of facility £ [0, 3d] 
where d is the number of clients. 


NOTE : Time is noted in CPU clocks, where clocks per second is equal to 
1000000. 
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Result Tables 1 


f 

d 

Primal-dual 

Linear time 
heuristic 

3 

10 

31 

69 

3 

50 

288 

261 

3 

100 

962 

472 

3 

500 

13587 

2217 

3 

1000 

26332 

4820 


Toal Time vs Demand (F=3) 
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Figure 56: Total time vs Demand (F = 3) 
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6000 



Figure 57: Clocks for linear time heuristic vs Demand (f=3) 
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Result Tables 2 


f 

d 

Primal-dual 

Linear time 
heuristic 

5 

50 

532 

1060 

5 

100 

1508 

1994 

5 

500 

15351 

8880 

5 

1000 

48834 

10883 

5 

2000 

187849 

19253 
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Figure 58: Total time vs Demand (F = 5) 
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Figure 59: Clocks for linear time heuristic vs Demand (f=5) 
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Result Tables 3 


f 

d 

Primal-dual 

Linear time 
heuristic 

8 

50 

618 

7385 

8 

100 

2199 

14113 

8 

500 

27781 

31999 

8 

1000 

79929 

56213 

8 

2000 

301885 

106905 


Toal Time vs Demand (F=8) 
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Figure 60: Total time vs Demand (F = 8) 
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Figure 61: Clocks for linear time heuristic vs Demand (f=8) 
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Result Tables 4 


f 

d 

Primal-dual 

Linear time 
heuristic 

10 

50 

735 

19704 

10 

100 

2451 

31716 

10 

500 

31085 

114011 

10 

1000 

106836 

235199 

10 

2000 

350755 

438162 

10 

5000 

1998171 

1057722 


Toal Time vs Demand (F=10) 
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Figure 62: Total time vs Demand (F = 10) 
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Figure 63: Clocks for linear time heuristic vs Demand (f=10) 
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Result Tables 5 


f 

d 

Primal-dual 

Linear time 
heuristic 

13 

50 

921 

117127 

13 

100 

2205 

204473 

13 

500 

33767 

917797 

13 

1000 

125070 

1808675 

13 

3000 

527109 

3603387 

13 

6000 

3773695 

10775159 


Toal Time vs Demand (F=13) 
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Figure 64: Total time vs Demand (F = 13) 
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Figure 65: Clocks for linear time heuristic vs Demand (f=13) 
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Result Tables 6 


f 

d 

Primal-dual 

Linear time 
heuristic 

15 

100 

3961 

838398 

15 

250 

14158 

1982816 

15 

500 

40603 

3859993 

15 

1000 

157742 

7688321 

15 

3000 

1278692 

22691206 

15 

6000 

4587448 

45293030 


Toal Time vs Demand (F=15) 
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Figure 66: Total time vs Demand (F = 15) 
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Figure 67: Clocks for linear time heuristic vs Demand (f=15) 
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Result Tables 7 


f 

d 

Primal-dual 

Linear time 
heuristic 

17 

100 

4201 

3490286 

17 

500 

45072 

16227801 

17 

1000 

183113 

32105374 

17 

3000 

1562383 

103233142 

17 

6000 

5294236 

205140190 


Toal Time vs Demand (F=17) 
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Figure 68: Total time vs Demand (F = 17) 
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Figure 69: Clocks for linear time heuristic vs Demand (f=17) 
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Result Tables 8 


f 

d 

Primal-dual 

Linear time 
heuristic 

20 

100 

2181 

30593379 

20 

500 

58569 

138222841 

20 

1000 

225908 

NIL 

20 

3000 

1912087 

NIL 

20 

6000 

6102226 

NIL 


Toal Time vs Demand (!F=20) 
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Figure 70: Total time vs Demand (F = 20) 


Analysis 

For / = 3 and / = 5 the total time taken by the primal dual scheme is 
higher for most of the outputs, for / = 8 and / = 10 a clear intersection is 
visible, for demand lower than that, total time taken by primal dual approach 
is lower and for demands higher than that, time taken is higher. As we in¬ 
crease the number of facilities, the time taken by the linear time heuristic 
is higher than that of primal dual approach for all the demand values, and 
the difference keeps growing as the number of facilities increase. This is so 
because the time complexity of the primal dual approach is polynomial with 
respect to both, demand d and number of facilities /, whereas it is linear with 
respect to demand d and exponential with respect to number of facilities / 
for linear time heuristic. Hence as the number of facilities increase, the com¬ 
putations required by the linear time approach overtakes that of primal dual 
approach and although gives a optimal solution, it becomes computationally 
too time consuming for practical use. 
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Analyzing the graphs for total time vs demand for linear time heuristic, it 
can be easily observed that total computational time required by the heuristic 
is linearly dependent on the demand size d. Which certainly is the case as 
the time complexity for this heuristic is 0(d) for constant value of /, where 
d is the number of clients. 

The graph given below is the measure of clocks required to compute the so¬ 
lution by the linear time heuristic for demand (d) = 100 for different number 
of facilities. By analyzing the graph, it can be easily observed that the com¬ 
putational time is exponentially proportional to number of facilities. Which 
can be seen in the time complexity of the heuristic that is 0(df 2-f -1 ) 


Total time vs number of Facilities (d=100) 
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Figure 71: Total time vs number of facilities (d = 100) 


Figure 72 shows the total time in clocks vs Demand graph for different 
number of facilities. One can observe how huge the change in slope is with 
just a slight increase in the value of /. Whereas individually, every curve 
is a straight line signifying the linear dependence of time with demand and 
massive increase in slope is the result of exponential dependence of time with 
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number of facilities /. And so much that lines for/ = 3,/ = 5,/ = 8 are 
not even visible. 
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Figure 72: Total time vs Demand for different number of facilities 


Conclusion: For random distribution of clients and facilities in a plane, a 
better solution can be obtained by using primal dual algorithm than using 
nearest neighbour approach if the number of facilities is large, more than 
250. Where as for small number of facilities nearest neighbour approach 
works better for large number of clients, but primal dual approach works 
better if the number of clients is also kept small. 

It is practically better to use the linear time heuristic over the primal-dual 
method provided the number of facilities is as low as 10 — 12, getting optimal 
results, instead of an approximate answer. For more than 13 facilities, is 
becomes computationally challenging and inefficient to compute an optimal 
solution by the linear time algorithm, preferring the primal-dual approach. 
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